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We evaluate baryon wave functions in both the Coulomb and Landau gauge in lattice QCD. These 
are constructed from quark propagators calculated with the overlap Dirac operator on quenched 
gauge configurations at /3 = 6. By comparing baryon states that differ in their diquark content, we 
find evidence for enhanced correlation in the scalar diquark channel, as favored by quark models. 
We also summarize earlier results for diquark masses in the Landau gauge, casting them in a form 
more easily compared with subsequent studies. 
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I. INTRODUCTION 

The notion of a diquark is nearly as old as that of 
quarks themselves and has been invoked to explain many 
aspects of hadron phenomenology (see [l[ for a review). 
Most generally, a diquark is any two quark system, but 
the term is more often taken to denote two correlated 
quarks in a particular representation of flavor and spin. 
In QCD-inspired quark models the color-hyperfine 
interaction gives rise to attraction in the spin singlet, 
SU(3)-flavor anti-triplet channel, a configuration known 
as a scalar diquark or more evocatively as a "good" di- 
quark. In contrast, the spin triplet, flavor sextet channel 
is repulsive, and the associated axial vector or "bad" di- 
quark is disfavored. Note that in this discussion and the 
rest of the paper, we only consider positive-parity di- 
quarks in the 3 of color, as would describe two valence 
quarks in a baryon. While one may write down diquark 
operators symmetric in color, all evidence points toward 
their being energetically disfavored. 

In recent years, diquarks have received increased atten- 
tion in light of the possible existence of exotic states such 
as the + , as diquark models make definite predictions 
for their properties Q. The status of the + remains 
uncertain (see [3| for a recent review of the experimental 
situation), but it serves to remind one of the relative lack 
of other exotics naively allowed by QCD, a scarcity that 
may largely be explained if diquark correlations play an 
important role in hadron structure Q. 

Ideally, issues such as these should be addressed by 
direct appeal to the fundamental theory. The lattice is 
the principal calculational framework for nonperturba- 
tive QCD and has been brought to bear on the question 
of diquarks in several recent studies (we set aside direct 
searches for exotic states). Perhaps the most straightfor- 
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ward approach is to construct a diquark two-point func- 
tion and consider its fall-off in time, as one does to ex- 
tract hadron masses. A diquark by itself is not a color 
singlet, however, and so one must either fix the gauge 
or introduce an additional source of color. The former 
approach was first pursued in [|| , where diquark correla- 
tors were calculated with Wilson fermions in the Landau 
gauge. More recently, we presented a similar investiga- 
tion in Q with overlap fermions at significantly lighter 
quark masses. By comparing the effective mass of the 
diquark with that of its constituent quarks, the scalar 
diquark was found to be bound in the limit of vanishing 
quark mass. In Section IIVI below, we briefly summa- 
rize these results in order to give values for mass split- 
tings that may be more easily compared with subsequent 
studies. In the second approach, one constructs a gauge 
invariant object by contracting the free color index of 
the diquark at source and sink with a Wilson line, serv- 
ing as a static quark [H, [§] . This allows one to extract 
diquark mass differences, in qualitative agreement with 
the fixed-gauge approach. See also (iq |. where point-to- 
point baryon correlators containing various diquarks are 
compared to those in the free theory. 

While useful, such mass determinations provide lim- 
ited information about the nature of diquark correlations. 
In this work, we directly investigate spatial correlations 
among quarks in baryons by calculating baryon wave 
functions on the lattice. At least two natural formalisms 
exist for defining what is meant by a "wave function." 
The one pursued here begins with a standard baryon cor- 
relator and involves displacing quarks at the sink. This 
function of quark displacements is then evaluated in a 
fixed gauge [11| • A very early study of such wave func- 
tions ma y b e found in [121 ] and more complete investiga- 
tions in jl3l . 03]. These treat only a subset of all pos- 
sible quark displacements and are largely motivated by 
a desire for improved interpolating operators for spec- 
troscopy. Nevertheless, and although not emphasized, 
the nucleon wave function parametrized in [13l . |14| ] does 
exhibit characteristics attributable to diquark effects, in 
particular a negative charge radius for the neutron. 
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An alternative definition of a hadronic wave func- 
tion is that provided by the density-density correla- 
tor method [TBI, [lfl [l7]|- A recent addition to the 
body of work treating such correlation functions for 
baryons [HI, 0, [2(| may be found in [2lJ , where the fo- 
cus is on possible deformations arising from spin-orbit 
coupling. Finally, a very recent study 0] employs the 
density-density correlator technique to examine the wave 
function of a diquark constrained to a spherical shell 
about a static quark, a gauge-invariant setup mentioned 
above in the context of diquark mass differences. By 
fitting to an exponential ansatz, the authors of Q find 
a large, but finite, radius for the scalar diquark in this 
environment. 

In this work, we present the first detailed study of di- 
quark correlations in physical baryons (with all quark 
masses finite). We consider all possible displacements of 
the three quarks and calculate wave functions in both 
the Coulomb and Landau gauges. By directly compar- 
ing wave functions of disparate states and calculating 
ratios of mean quark separations, we find evidence of 
enhanced correlation in the scalar diquark channel. We 
work in quenched QCD and employ the overlap Dirac op- 
erator [H, [H, HH El| in our calculation, a discretization 
which preserves chiral symmetry on the lattice [H, [53] 
and is thereby closest to the continuum formulation. 

The paper is organized as follows. In Section [Til we 
provide details of our calculation and describe the corre- 
lation functions and states that we study. In Section ITTTI 
and its subsections, we present and compare our baryon 
wave functions and from them calculate values for mean 
quark separations. Finally, in Section HVl we calculate 
diquark mass differences from data first presented in [?J . 

II. DETAILS OF THE CALCULATION 

This study is one in a series employing the overlap 
Dirac operator on a large lattice. Results for meson and 
baryon spectra, as well as meson wave functions, diquark 
correlators, and other observables were presented in [7j. 
In [281 ]. we calculated matrix elements relevant for kaon 
physics in the standard model and beyond with a care- 
ful treatment of nonperturbative renormalization in the 
RI/MOM scheme. We direct the reader to [7| for a discus- 
sion of the many advantages of the overlap discretization 
as well as for details of our implementation beyond those 
given here. 

The overlap Dirac operator describing a massless 
quark [23[ is given by 

a V \fxJx ) 

where X = Dw — p/a is the Wilson Dirac operator with 
mass —p/a. It follows that inversion of the overlap oper- 
ator requires the repeated calculation of 1/ yXJX. This 
is accomplished with polynomial or rational function ap- 
proximations and is very demanding computationally. 



An unquenched calculation on a lattice as large as ours 
would be beyond the capability of presently available re- 
sources. We therefore work in the quenched approxima- 
tion and note that prior experience with Wilson fcrmions 
has shown hadronic wave functions of the type we study 
to be largely unaffected by quenching (l4| . 

We employ the Wilson gauge action with (3 = 6 on a 
lattice of size 18 3 x 64. This gives an inverse lattice spac- 
ing a -1 of 2.12 GeV [H, [3(| on the basis of the Sornmer 
scale defined by 7oF(r ) = 1.65 with r = 0.5 fm [3l| . 
One hundred independent gauge configurations were gen- 
erated and then fixed to the Landau gauge before invert- 
ing the Dirac operator. The negative mass parameter in 
the definition of the overlap was set to p = 1.4 in order 
to maximize locality [32| . Quark propagators were calcu- 
lated from a point source for all color-spin combinations 
with a conjugate gradient multimass solver for bare quark 
masses am q = 0.03,0.04,0.06,0.08,0.10,0.25,0.50,0.75. 
For reference, corresponding values of the pion mass as 
well as baryon masses are given in Table UJ [7| . 

We now describe the states that we study. We work 
in a Dirac basis where 74 is diagonal and utilize "non- 
relativistic" wave functions involving only either upper 
or lower spinor components. Labeling the three quarks 
u, d, s for convenience, we give the spin structure of the 
states of interest in Table UTI in a transparent notation. In 
the SU(3) classification, these correspond to the octet A 
and E states and the decuplet E*. In Section [TTTl we will 
find it most illuminating to compare the A to the E*; in 
the former, the u and d are in a spin-0 "good diquark" 
configuration, while in the latter they are in a spin-1. 
The octet E is a cousin of the nucleon in which the pairs 
u, s and d, s are in superpositions of spin-0 and spin-1. 

For a given state, with spin structure as given in the 
table, we construct a zero-momentum correlator 

G(f u ,f d ,t) = y^(u(r + r u ,t)d(f+ f d ,t)s(f,t) 

r 

x u(0, 0)J(0, 0)s(0, 0)) . (2) 

Here color indices are implicit and are contracted with 
the antisymmetric tensor at source and sink. We com- 
bine correlators for the (two for the octet, four for the 
decuplet) spin states distinguished by J z . Finally, to 
the forward-propagating correlators constructed with up- 
per spinor components we add correlators propagating in 
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aMio 


0.03 


0.219(3) 


0.63(2) 


0.75(3) 


0.04 


0.247(2) 


0.66(2) 


0.78(2) 


0.06 


0.297(2) 


0.714(11) 


0.82(2) 


0.08 


0.340(2) 


0.763(9) 


0.868(12) 


0.10 


0.3803(14) 


0.810(7) 


0.909(10) 



TABLE I: Masses, in lattice units, of the lightest pseudo-scalar 
meson, octet baryon (e.g. nucleon), and decuplet baryon for 
quarks of equal mass am q . Quoted errors are statistical only. 
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TABLE II: Baryon states. 



the backward time direction that have been constructed 
with lower components. We thereby double our statis- 
tics while ensuring that only the desired positive-parity 
states are excited from the vacuum. 

Since the quarks at the sink may be taken to be 
at distinct spatial sites, Eq. @ is only well-defined if 
we specify the gauge. In Section IIIII we evaluate this 
correlation function in both the Coulomb and Landau 
gauges. Coulomb gauge-fixing was performed using sim- 
ulated annealing, starting from gauge configurations al- 
ready fixed to the Landau gauge. At sufficiently large 
times, G(f u ,fd,t) settles into a spatial profile that is in- 
dependent of t up to normalization. We refer to this 
profile as the "wave function," 



*(r„,rd) 



G(f u ,fd,t) 



(3) 



This zero-momentum wave function in general depends 
on two 3-vectors, i.e. six numbers. As discussed in the 
next section, however, we only resolve a dependence on 
separations between pairs of quarks, and it is therefore 
effectively a function of a triangle, parametrized by three 
numbers. For the purpose of displaying the wave func- 
tion, we adopt the geometry shown in Fig. [TJ Here z is 
the distance between the quark labeled by s and the cen- 
ter of mass of u and d. The axis determined by z is taken 
to establish a coordinate system in which we specify the 




position (x, y) of u with respect to the center of mass. We 
note that the states we consider are all symmetric under 
interchange of the positions of u and d. In summary, the 
wave function in these coordinates is given by 
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(4) 



where we have normalized the amplitude to unity where 
all three quarks are at the same site and have defined 
the delta function on the lattice taking into account the 
multiplicity of the sites. In constructing the wave func- 
tion, we only consider configurations of the quarks where 
no two are separated by more than half the length of the 
lattice (L/2 = 9a). 

We conclude this section with some final details of our 
implementation. First, we note that it might be advan- 
tageous to replace the point source in Eq. ([2]) with an 
extended operator that better overlaps the desired state. 
We were constrained in our calculation, however, by the 
fact that point-source propagators were required for stud- 
ies of nonperturbative renormalization and weak matrix 
elements; the calculation of an additional set of smeared- 
source propagators was deemed too costly to be worth- 
while. 

The wave functions we present in Section [IIII were cal- 
culated by summing over all possible positions of the 
three quarks. Since each sum is over 18 3 sites, this in- 
volves a nontrivial amount of work. We were able to 
greatly speed up the calculation, however, by employ- 
ing a fast Fourier transform and utilizing the convolution 
theorem to eliminate one of the summations. A related 
issue is the large amount of data that would have to be 
stored to capture all possible relative displacements of the 
quarks (i.e. all possible embeddings of a triangle in the 
lattice). This was avoided by adopting the parametriza- 
tion described above and building a histogram in the 
x, y, z coordinates with linear interpolation. The bin size 
was taken to be 0.225a in x, y and 0.45a in z, sufficiently 
small that the mean quark separations presented in Sec- 
tion IIII Bl arc unbiased, as confirmed by examining the 
totally symmetric E* state. 

Computations were performed with shared memory 
code on IBM p690 systems at Boston University and 
NCSA. 



III. BARYON WAVE FUNCTIONS 



A. Wave functions 



FIG. 1: Geometry for visualizing the wave functions. 



Visualization of the wave functions will prove to be 
quite useful for discerning differences in spatial correla- 



FIG. 2: (Color online) Wave function of the A evaluated at 
t — 10a in the Landau gauge, for z — 2.25a. 



tions between states. As a first step, we must choose the 
time t at which to evaluate the wave function. For small 
times, the correlator G(r u ,rd,t) in Eq. @ is dominated 
by excited states. It is therefore necessary to take t suffi- 
ciently large that the spatial profile has settled into that 
of the ground state. We find that for the states we study, 
the wave function has settled by t — 8a, in agreement 
with what was observed for effective masses when calcu- 
lating baryon spectra Q. We conservatively take t = 10a 
in the remainder of this paper. 

For plotting purposes, we fix z, the distance between 
the center of mass of the first two quarks and the posi- 
tion of the third. In Fig. [21 we plot the A wave func- 
tion in the Landau gauge as a function of x, y for one 
such z separation. All three quark masses are taken to 
be am q — 0.03, the lightest available value. The corre- 
sponding wave function in the Coulomb gauge is plotted 
in Fig. [3J Recall that the wave functions have been nor- 
malized to 1 where all three quarks are at the same site 
(x = y = z = 0, not shown). Figures [2] and [3J exhibit 
the general property that Coulomb-gauge wave functions 
are less broad and better contained in the lattice volume 
than those calculated in the Landau gauge, in agreement 
with [l3|, [3]. We will focus on Coulomb-gauge wave 
functions in the remainder of this section. 

Statistical errors, which for clarity are not shown in 
Figs. are on the order of 6 to 10 percent (see Fig. [S] 
below). It is noteworthy that the overall amplitude of 
the wave function tends to vary configuration by con- 
figuration while it maintains the same basic shape. In 
other words, if the wave function at its peak (always at 
r u = = 0) is found to be larger than average on a given 
gauge configuration, it is likely to be larger at all other 
quark displacements on that configuration. In (l3l. [T^|. 
this effect was taken as motivation to normalize the wave 



FIG. 3: (Color online) Wave function of the A evaluated at 
t — 10a in the Coulomb gauge, for z = 2.25a. 



functions on a per-configuration basis. While effective, 
this approach is difficult to justify from a field-theoretic 
perspective and we do not pursue it here. We note, how- 
ever, that when comparing the properties of various wave 
functions quantitatively, such contributions to the errors 
often cancel, as we find for mean quark separations in 
the next section. 

By parametrizing our wave functions in terms of rela- 
tive separations, without regard to orientation, we have 
implicitly assumed isotropy. Of course, one recognizes 
that there is a preferred direction, the z-direction of the 
lattice (not to be confused with our z coordinate) with re- 
spect to which the ^-component of spin is defined. To test 
for the possible presence of spin-orbit coupling, we added 
a fourth dimension to our histogram with the new vari- 
able being the projection of the vector whose length we 
call "z" along the z-direction of the lattice. We then con- 
structed a wave function with definite J z and looked for 
a dependence on this variable. Within errors, we found 
no evidence for such a dependence. We conclude that 
the effects of spin-orbit coupling, if present, are below 
the statistical limits of our calculation. 

We come now to the main point of interest. For the 
A state whose wave function is plotted in Fig. [31 the it, d 
quarks are in the spin-0, "good diquark" configuration. 
We would like to compare this wave function to that of 
the £*, where the two quarks are in the spin-1 configura- 
tion. In Fig. [Hand Fig.0 we plot the two wave functions 
together for two different z separations. A cross-section 
of Fig. Q] with y — is shown in Fig. [5] to give an indi- 
cation of the errors. As predicted in the literature, we 
note significantly stronger clustering when u, d are in the 
good diquark configuration. This feature is independent 
of z. 

In this section, we have presented results for baryons 
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FIG. 4: (Color online) Comparison of A (red) and S* 
(broader, in green) wave functions in the Coulomb gauge, for 
2 = 2.25a. 



FIG. 6: (Color online) Cross-sections of A and E* wave func- 
tions with z — 2.25a and y = 0. 
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two light particles bound by two-body interactions. 



B. Mean quark separations 

From our wave functions, we calculate the mean square 
separation between the u and d quarks in the natural 
way: 



(k-^heEi*^ 



rd)\ \r u -r d \ 



(5) 



Similarly, noting that our coordinates are defined such 
that r s — 0, 



(6) 



FIG. 5: (Color online) Comparison of A (red) and E* 
(broader, in green) wave functions in the Coulomb gauge, for 
z = 4.50a. 



where the three quarks are taken to be degenerate in 
mass with am q = 0.03. fn the next section, we will eval- 
uate the effect of increasing this mass. One may also 
consider baryons where the s quark is taken to be signif- 
icantly heavier than the others. None of the qualitative 
features are changed, but we do observe a slight tendency 
for the mean separation between the two light quarks to 
be larger than that between one of the light quarks and 
the heavy quark, when considering the otherwise sym- 
metric S* state. This is a purely kinematic effect that 
would apply even in a classical system of one heavy and 



Since our baryons reside in a finite volume, such mean 
separations must be interpreted with care, fn particu- 
lar, we do not take into account the tails of the wave 
functions that extend into adjacent cells of our periodic 
lattice, nor do we remove those that impinge from them. 
To do so would require that we model and fit the numer- 
ical wave functions, fn contrast, the separations that we 
calculate follow directly from the data, fn the large vol- 
ume limit, these separations would converge to definite 
values. In our finite volume, they provide a rough quan- 
titative estimate of the clustering observed in the scalar 
diquark channel and of the dependence of such clustering 
on quark mass. To the extent that finite volume effects 
are present, they are expected only to weaken correla- 
tions. In Tables Ifffl and IfVI in the appendix, we collect 
root mean square separations for the various states in the 
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FIG. 7: (Color online) Ratios of RMS quark separations as a 
function of bare quark mass am q . 
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FIG. 8: (Color online) Constituent quark and diquark masses 
in the Landau gauge, taken from [7j. 



two gauges. Errors have been calculated via the boot- 
strap method with 500 samples. 

In QCD-inspired quark models 0, the Hamiltonian 
generally includes a term of the form 

H = o s cV — - — Si-sj (7) 

i<j 

that is attractive in the spin singlet channel. Here Sj is 
the spin and rrii the (constituent) mass of the ith quark, 
and c is a constant. It follows that the strength of the 
interaction increases as quark masses decrease. For the 
purpose of quantifying the mass dependence of our wave 
functions, we define a ratio of RMS separations, 




As noted earlier, fixed-gauge wave functions are gener- 
ally broader in the Landau gauge than in the Coulomb 
gauge. For example, for the A state with am q — 0.03 
in the Landau gauge, we find \J (\f u — rd\ 2 ) = 5.63(6), 
as compared to 5.17(9) in the Coulomb gauge. Remark- 
ably, however, ratios of separations appear to be rather 
independent of gauge. 

In Fig. [71 we plot such ratios in both gauges for the 
two octet states and all available quark masses. By con- 
struction, the decuplet E* is totally symmetric; the corre- 
sponding ratio is exactly one and would lie on the dotted 
line in the figure. We recall that in the A, the u, d are 
in the spin-0 configuration while in the E, the u, s and 
u,d are in superpositions of spin-0 and spin-1. Again, 
errors have been calculated with the bootstrap and leave 
little doubt that spatial correlations are enhanced in the 
scalar diquark channel. We also observe that the effect 
strengthens markedly at the lightest masses. 



IV. DIQUARK MASSES 

In the above sections, we have observed diquark ef- 
fects via spatial correlations in baryon wave functions. 
An alternative approach for investigating diquarks on the 
lattice is to construct diquark-diquark correlators and fit 
their decay in Euclidean time in terms of an effective "di- 
quark mass" [6j]. This is not a gauge-invariant concept 
and such a parameter cannot be interpreted as the mass 
of a physical state, but it may nevertheless give some in- 
dication of the relative strength of binding. In we 
presented results for diquark masses calculated in the 
Landau gauge from correlators of the form 

G M = 2<etf*t«i(»*.*)d*(f,i)e«'A'«i'(ff,0)(ijk/(ff,0)) , 

r 

(9) 

where the indices label color, and implicit spin indices are 
assigned such that u,d are in either the spin-0 or spin- 
1 configuration. In Fig. [5J we reproduce a plot taken 
from [7J, showing the dependence of diquark masses on 
quark mass. Also included is the "constituent quark 
mass," determined by performing a fit to the quark prop- 
agator in the Landau gauge. If one takes seriously this 
"constituent mass" interpretation, it appears that the 
scalar diquark may be bound in the limit of vanishing 
quark mass. Here we expand on these earlier results in 
two ways, by utilizing non-local sinks in the construction 
of the correlators and by reporting values for mass split- 
tings, with errors taking into account correlations in the 
data. 

In Q, diquark masses were extracted from point-to- 
point correlators. As discussed in Section UH we remain 
constrained to point sources, but we are free to use ex- 
tended sinks, following the approach applied to mesons 
in Q. A natural choice for this purpose is to use the 
diquark analogue of the baryon wave functions presented 
earlier, again calculated in the Landau gauge. In par- 
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FIG. 9: (Color online) Effective mass of the scalar diquark as 
a function of time, calculated with both point and extended 
sinks. The effective quark mass is also shown. 
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FIG. 10: (Color online) Effective mass of the vector diquark 
as a function of time, calculated with both point and extended 
sinks. The effective quark mass is also shown. 



ticular, we generalize Eq. ([9]) to allow for a separation r 
between quarks at the sink, 

G(r,t) = (eijkUj(f u ,t)d k (r d , t)e i j>i t >'Uj>0, 0) 

Pu,r d 

xd k ,(6,0))5(r-\r u -r d \). (10) 

We take this function of r at t = 10a to define the wave 
function, (p(r) — G(r, 10a), calculating such a y(r) for 
each state and quark mass of interest. Finally, from these 
we construct an extended-sink correlator, 



G cx t(i) =5>(r)G(r,i), 



(11) 



whose fall-off yields the desired diquark mass. 

In Fig. [9l we plot the effective mass of the scalar 
diquark as a function of time, given by aM c g = 
—ln[G(i)/G(t — a)], for both point-sink and extended- 
sink correlators. The bare quark mass is am q = 0.03, 
our lightest value. The corresponding plot for the vector 
diquark is shown in Fig. [TO] We find that both point and 
extended sink correlators display the same asymptotic ef- 
fective mass, giving us confidence that the observed rate 
of exponential decay may be interpreted as the "mass" of 
the corresponding diquark in the Landau gauge. Results 
at heavier quark masses exhibit similar behavior. For the 
results that follow, we fit the extended-sink correlators in 
the region 11 < t/a < 14 and calculate statistical errors 
by bootstrap. 

We first consider the "binding" of the scalar diquark 
with respect to the combined mass of two "constituent 
quarks," indicated earlier in In Fig. [TTJ we plot the 
difference between the scalar diquark mass Ms=o and 
twice the constituent quark mass Mq. A naive linear 
extrapolation gives a(Ms=o — 2Mq) = —0.10(4) in the 
chiral limit. 

A more robust feature than the binding of the scalar 
diquark is the large splitting between it and the vector 
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FIG. 11: Difference between the mass of the scalar diquark 
and twice the constituent quark mass, as a function of bare 
quark mass. 



diquark. We plot our results for this mass difference in 
Fig. Q21 A linear extrapolation to the chiral limit gives 
a{M s =i-M s =o) = 0.077(35). Taking a = 2.12 GeV from 
the Sommer scale, we find M s =i-M s =o = 162(75) MeV, 
where the error is statistical only. This splitting has also 
been calculated in HQ in a gauge- invariant setup where 
the free color index of the diquark operator at source and 
sink is contracted with a Wilson line. Equivalently, this 
scheme corresponds to evaluating the diquark correlator 
in a temporal gauge, in which the temporal Wilson line 
reduces to the identity. Our different choice of gauge 
does not allow a direct comparison, but we note that 
the splitting is universally found to be positive and that 
it might be interesting to further investigate the gauge 
dependence of this quantity. 



0. 1 r 



. 0.05 ■ 



0.02 0.04 0.06 0.( 

am 

q 



0.1 



FIG. 12: Mass splitting between the scalar and vector di- 
quark, as a function of bare quark mass. 



V. CONCLUSIONS 

In this work, we evaluated baryon wave functions in 
the Coulomb and Landau gauges and compared them on 
the basis of their diquark content. We found that spatial 
correlations were significantly enhanced between quarks 
in the scalar diquark configuration as compared to the 
vector diquark. Finally, we presented results for effective 
mass differences between diquark states calculated in the 
Landau gauge. 

We acknowledge that our calculation suffers from lim- 
itations of the quenched approximation and the manifest 
gauge-dependence of our wave functions. It is encour- 
aging, however, that enhanced correlations were equally 
pronounced in both gauges. It is also noteworthy that 
in all cases, diquark effects were found to become more 
pronounced as quark masses were decreased. This mo- 
tivates further, preferably unquenched, calculations at 
lighter masses. 
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APPENDIX: TABLES OF QUARK 
SEPARATIONS 

In Tables IIIII and IIV1 we provide the root mean square 
separation between quarks in the various states, calcu- 
lated as described in Sec. IIII Bl 
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TABLE III: RMS separation yj (\n - rj\ 2 )/a, in lattice units, 
between quarks of flavor i and j as a function of bare quark 
mass, from baryon wave functions evaluated at t = 10a in the 
Landau gauge. 
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TABLE IV: RMS separation ^/(|fi - r j | 2 )/a, in lattice units, 
between quarks of flavor i and j as a function of bare quark 
mass, from baryon wave functions evaluated at t = 10a in the 
Coulomb gauge. 
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